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ABSTRACT 

This paper presents long slit spectral maps of the bi-cone shaped extended narrow 
line region (ENLR) in the Seyfert galaxy NGC 2992. We investigate the physical prop- 
erties of the ENLR via emission line diagnostics, and compare the observations to shock 
and photoionization models for the excitation mechanism of the gas. The line ratios 
vary as a function of position in the ENLR, and the loci of the observed points on line 
ratio diagrams are shown to be most consistent with shock+precursor model grids. We 
consider the energetics of a nuclear ionizing source for the ENLR, and perform the g-test 
in which the rate of ionizing photons from the nucleus is inferred from measurements of 
the density and ionization parameter. The g-test is shown to be invalid in the case of 
NGC 2992 because of the limitations of the [S II] A6717/A6731 density diagnostic. The 
excitation of the gas is shown to be broadly consistent with the kinematics, with higher 
[N II] A6583/Ha present in the more dynamically active region. We also show that the 
pressure associated with the X-ray emitting plasma may provide a large fraction of the 
pressure required to power the ENLR via shocks. 



Subject headings: galaxies: active — galaxies: individual (NGC 2992) — galaxies: 
Seyfert 



1. INTRODUCTION 

The narrow line regions (NLR) of active galaxies are the largest structures which directly 
derive their energy input from the active galactic nucleus (AGN). In many nearby Seyfert galaxies 
the NLR are spatially resolved, and often extend over distances of several kpc, so that the energy 
input processes may be studied in detail. Such spatially-resolved NLR are often referred to as 
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extended narrow line regions (ENLR) flUnger et al.| 1987). In both Seyfert 1 and Seyfert 2 galaxies 
the ENLR appears as a filamentary system of high excitation gas and in Seyfert 2 galaxies this gas 
is often confined to a bi-cone structure centred on the nucleus, the so called "ionization cones" . In 
other cases, the morphology is more complex, and bent jet-like structures are sometimes observed. 



According to the model developed some twenty years ago ( [Koski 1978; Ferland & Netzer 1983; 
Stasinska] 1984), and which remains the dominant paradigm ( |Robinson et al. 1987, [Viegas fc dej 
Gouveia Dal Pino] 1992), the NLR are excited by photoionization by a powerful (approximately 
power-law) non-thermal UV continuum produced by the AGN. This scenario appeals to the unified 
schemes for AGN in which the nucleus and the broad-line region are surrounded by an optically 
thick torus ( |Antonucci fc Miller| 1985). In this picture the differences between Seyfert types 1 and 
2 are ascribed to an orientation effect in which type 1 Seyferts are oriented so that the line of sight 
falls within the opening angle of the cones, while in type 2 Seyferts the line of sight to the nuclear 
region is obscured by the torus. This scenario explains the conical ENLR observed in some type 2 
Seyferts as due to shadowing of the ionizing radiation from the nucleus by the absorbing torus. 

Shocks may also be important for the morphology, kinematics and excitation of the NLR in 
Seyfert galaxies. In some Seyferts the morphology of the emission line gas is directly related to 
that of the radio emission, providing evidence for dynamical interaction between the radio plasma 
and the ambient interstellar medium of the host galaxy. For example, high resolution narrow band 
imaging with HST has revealed bow-shock shaped emission line regions around radio lobes, and 
jet-like emission line regions in sources with jet-like radio structures ( |Axon et al. 1998, Pogge 1996, 
and references therein). 

Shocks probably dominate the physics of the large scale outflows observed in Seyfert galaxies. 



In a series of papers, ( Colbert et al. 1996a, Colbert et al. 1996b, Colbert et al. 1998) Colbert and 
collaborators show that optical emission suggestive of large scale outflows is present in ^ j of their 
sample of 22 edge-on Seyfert galaxies, ^ i of the sample have kiloparsec scale radio structures, 



and 3 of their sample, including NGC 2992, have similarly extended soft X-ray emission. Colbert 



et al. (1998) favour an interpretation in which these large scale outflows originate as AGN driven 
non-thermal jets that entrain the surrounding thermal material, and heat it by shocks over kpc 
spatial scales. An alternative possibility is that the X-ray gas results from starburst driven winds 
from the circumnuclear region. 

Seyfert galaxies are radio quiet, and simple equipartition arguments have long suggested that 
the pressure in the relativistic plasma responsible for the radio wavelength synchrotron emission 



is inadequate to drive high velocity shocks into the interstellar medium ( Wilson, Ward, &: Haniff 
1988). However, this picture may be wrong if mass entrainment into the jet occurs close enough 



to the source. Recently Bicknell et al. (1998) has shown that although the radio jets of Seyferts 
start off mildly relativistic, over most of their length the internal energies of Seyfert jets must 
be dominated by thermal plasma. The combined pressure of the thermal and relativistic gas is 
then sufficient for the jets to power the optical emission of the NLR. This idea finds observational 
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support from the work of polbert et al.| (1998) which suggests that the large scale outflows in 
Seyfert galaxies are predominantly winds of thermal X-ray emitting gas. This thermal phase will 



dominate the pressure, consistent with the requirements of the theory described in Bicknell et al. 
(1998). 

Models of shocks for Seyfert galaxies have considered various scenarios for the excitation of 
the gas. Bow-shock models developed by [Taylor, Dyson fc Axon| (1992) (which can account for the 
wide range of line profile shapes observed in the NLR) consider the excitation mechanism to be 
photoionization, whereby the shocks simply serve to compress the gas which is then illuminated 
by UV continuum from the nucleus. More recent models show that shocks themselves can provide 
an alternative mechanism for the ionization of the NLR based on the input of mechanical energy. 
These shock models advanced by [Sutherland, Bicknell &: Dopita| (1993) and developed in detail 
by Dopita fc Sutherland! (1996) show that high velocity shocks generate strong UV radiation in 



the high temperature cooling zone behind the shock front, the hardness of which depends on the 
shock velocity. This photon field is then available to produce a highly ionized precursor H II region, 
which emits a spectrum like that observed in the NLR. Dopita fe Sutherland! (1995) find that their 



shock+precursor models can account for the observed optical emission line ratios of a sample of 
Seyfert galaxies. This raised the possibility that such models may be generally applicable to narrow 
emission line regions. 

In order to find observational support in favour of either the shock or the photoionization 
model of the NLR, we have undertaken a detailed spectroscopic and dynamical study of the ENLR 
of a number of nearby Seyferts. Amongst these, the subject of this paper, NGC 2992 provides a 
particularly good test case to compare the predictions of the models with observations. NGC 2992 
is a gas-rich ( [Hutchmeier 1982; Sanders and Mirabel 1985) Sa type galaxy classified as having 



a 



Seyfert 1.9 active nucleus. It is inclined at ~70° to our line of sight and is crossed by a disturbed 
dust lane. It is interacting with NGC 2993 which is connected to the southern part of NGC 2992 
by a tidal tail with a projected length of 2.9'. The nuclear redshift is z=0.0077±0.0002 which gives 
angular scale 150 pc arcsec -1 assuming Ho=75 km s _1 Mpc -1 . 

It provides an interesting test case because it displays a striking and unusually extensive bi- 
conical ENLR which emerges almost perpendicularly from the plane of the galaxy. Narrow band 
images show the spatial extent of bright bi-cones is ~4 kpc, with the SE cone brighter than the 
NW cone. Fainter emission is extended up to 8 kpc. Figures la and lb show narrow band filter 
images obtained at the ESO New Technology Telescope (NTT). These are respectively [O III] A5007, 
and Ha+[N II] and are presented on a logarithmic stretch so as to bring out the fainter regions. 
Figure lc shows a blue continuum image of NGC 2992. 



The radio structure of NGC 2992 is described in detail in a paper by [Wehrle fc Morris] (1988) 



which includes a synopsis of previous radio observations. In summary NGC 2992 has 20 cm radio 
emission which extends approximately 25" along PA~160°, and has a one sided extension along 
PA~130°. At a wavelength of 6 cm there is visible a smaller scale figure-8-shaped structure, with 
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its long axis extending approximately 8" along PA~-26°, and centred on the nucleus. The figure-8 
may represent edge brightened bubbles of relativistic plasma or else magnetically-compressed arches, 
either of which are indicative of an expanding region of over-pressured relativistic plasma. However, 
comparison of the radio images with optical narrow band images shows no obvious correlation 
between the structures - see Figures Id, le and If where these radio maps are overlaid on our 
[O III] A5007, and continuum images. 

NGC 2992 was included in Durret and Bergeron's investigation of the physical properties of 



the ionized gas in emission line galaxies ( purret and Bergeron 1988, Durret 1990). In the earlier 



paper they present long slit observations of NGC 2992. The second paper describes the physical 
properties of the ionized gas and other general characteristics derived from their spectra. Here we 
briefly summarise their results. 

First, in the nuclear spectrum a broad component of Ha is detectable, but there is no corre- 
sponding broad H/3 component. NGC 2992 is therefore classified as an intermediate type Seyfert 



(1.9). This classification is confirmed from the shape of the far-IR to far-UV continuum (Boisson 



and Durret 1986). Second, Durret and Bergeron (1988) found the extranuclear spectra along PA = 
197° have similar line ratios to the nuclear region, and are characteristic of those found in Seyfert 2 
galaxies. This excludes the possibility that the ENLR is excited by OB stars. The excitation of 
the gas appears to decrease from the nucleus to the outer regions. The reddening of the inner 
regions is high, ranging up to ~3, and a similar reddening was estimated by |Malkan| (1983) 



from the [S II] AA(4069+4076)/(6717+6731) ratio. Third, using the Ha/H/? to correct line ratios 
for reddening purret and Bergeronj (1988) found a temperature of T e =19900K from the [O III] lines 
in the nuclear regions. 



The velocity field of NGC 2992 is complex. A kinematical study by Colina et al.| (1987) 



revealed blue asymmetries in the [O III] lines at the center of the galaxy, and showed that the 



nuclear region is dynamically decoupled from the off-nuclear regions. Colina et al. (1987) find 
some indication of radial outflow within a plane highly inclined to the galaxy plane, and that the 
kinematics of the SE region may be accounted for by tidal interaction with NGC 2993. Recently 



Marquez et al] (1998) have reported a kinematical study of NGC 2992 with the Ha and [O III] 



lines measured along nine position angles. They note that outflow of the gas within a conical 
envelope or on the surface of a hollow cone as one possible picture to account for the double peaked 
profiles and assymetries detected in several regions. However they use a simple kinematical model 
of circular rotation with a constant radial outflow, where the out-flowing velocities are considered 
as radial motions in a plane very close to that of the gas disk. While this model provides a good 
approximation to the observed kinematics we caution that such a model is not unique, and prefer 
to interpret the kinematics as a general outflow or wind, not necessarily constrained to a particular 



plane. This was reported briefly by us in Tsvetanov, Dopita & Allen (1995) and details will appear 



in a future paper. In summary, our higher spectral resolution study of the NGC 2992 NLR shows 
double component line profiles at most positions. One of the kinematical components follows a 
galactic rotation curve, while the other component is identified as a wind which is expanding out 
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of the plane of the galaxy with velocities up to 200 km s 1 . 

These kinematical studies generally assume the disk of NGC 2992 to be oriented with NW edge 
closest, in concordance with trailing spiral arms. In this orientation, with the axis of the bi-cone 
perpendicular to the disk, the SE cone is the closer and is directed towards us, consistent with the 
outflow interpretation. If on the other hand the spiral arms lead, then the SE edge of the galaxy 
and the NW cone would be closer. In this orientation the kinematics would indicate inflow instead 
of outflow. Other possibilities also exist if the bi-cone axis is highly inclined with respect to the 
galaxy disk. 



The X-ray spectral properties of NGC 2992 are described by Weaver et al. (1996) (and refer- 
ences therein) to have a significant absorbing column and a prominent Fe Ka fluorescence emission 
line. The X-ray flux varies by a factor of 2—3 on timescales of days to months. There is an excess 
of soft X-ray (0.1—2.0 keV) flux above that expected from the extrapolation of the 2—10 keV power 
law (r=1.7), and there is also some evidence for extended soft X-ray emission on a scale of ~2 
kpc. X-ray images of NGC 2992 taken with the ROSAT Position Sensitive Proportional Counter 



(PSPC), and High-Resolution Imager (HRI) are presented by |Colbert etaL| (1998). Soft X-ray (0.2 
— 2.4 keV) emission extends out to ~l' along PA~100°, and is roughly co-spatial with the 20 cm 
radio emission. 

Recently |Glass| (1997) have reported that the near infrared light-curve of NGC 2992 (monitored 
~3 times a year over the past 15 years, with an 9" aperture) shows an outburst in 1988. The 
amplitude of the event was greatest in the L(3.5 /im) band, which increased 102 mJy above the 
base level in less than 254 days, then decayed exponentially with an e-folding time of ~900 days. 
Glass| (1997) compare the total energy output of the event (6xl0 43 J) to the explosion energies 



of types I and II supernovae. They find that while a conventional un-obscured supernova does 
not provide a direct fit, the outburst may be explicable in terms of a model similar to that of 



Terlevich et al. (1992), which seeks to explain the characteristics of Seyfert nuclei in terms of 



multiple supernova explosions in dense H II regions. 

This paper presents the results of mapping the ENLR of NGC 2992 with long slit spectroscopy 
in order to investigate the detailed spatial dependence of line ratios and the physical properties 
with the objective of understanding the mode of excitation of the ENLR. The observations and 
data reductions are described in section 2, including the line flux measurements and representative 
spectra. Three distinct spectral regions of the ENLR are identified in section 3 in which we show 
the spatial variation of the line ratios. In section 4 we derive the physical properties of the emission 
line gas. The observations are compared to shock and photoionization model predictions with line 
ratio diagnostic diagrams in section 5. We discuss the g-test for photoionization from a central 
source in section 6. In section 7 we compare the excitation of the gas with the kinematics, and 
consider the possibility of a thermal pressure energy source for the excitation of the ENLR. 
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2. OBSERVATIONS AND DATA REDUCTION 

NGC 2992 was observed with the Double Beam Spectrograph (DBS) on the ANU 2.3m tele- 
scope at Siding Spring, Australia on March 5th and 9th 1994. Loral CCD detectors were used on 
both the red and blue arms of the spectrograph. The 158 line mm -1 grating was used in the red 
arm giving a dispersion of 4.1 A pixel -1 , a spectral resolution of 11.9 A and covering a wavelength 
range 5840-9725 A. The 300 line mm grating was used in the blue arm, providing a dispersion 
of 2.2 A pixel -1 , and yielding a 4.9 A spectral resolution over the wavelength range 3440—5570 
A. The slit width was fixed at 2" for all observations of NGC 2992, but was increased to 5" for 
observations of standard stars so as to ensure that slit losses were negligible. The spatial scale 
along the slit was the same in both red and blue arms; / /89 pixel -1 . 

Spectra were obtained for a total of nine slit positions all orientated at PA 30°. These were 
placed ~2" apart so as to efficiently cover the ENLR and provide a set of spatially independent 
spectra. A total of three 1000 s exposures were taken at each slit position. The FWHM of the 
standard star profiles along the wide slit provides an estimate of the seeing plus the guiding errors. 
At the beginning of the first night the median seeing was ~2" but improved to ~1'.'8 by the end of 
the night. The seeing in the middle of the second night was ~2'.'3. 

Data reduction was performed using the standard IRAF routines. Pixel to pixel variations 
in the CCD response were removed by appropriately normalised spectral flat fields derived from 
Quartz-Iodine lamp spectra, while a correction for the slit function was made by a low order fit to the 
illumination of the slit by the twilight sky. The DBS has very low vignetting, so this correction is not 
very important. The wavelength scale was established by fitting 5th order Chebychev polynomials 
to lines identified in Helium- Argon comparison lamp spectra obtained throughout the observations 
on both nights. 

Long-slit observations are usually affected by misalignments and optical distortions which cause 
the dispersion direction not to be exactly aligned along the lines of the CCD. Spectra of standard 
stars through a narrow slit revealed this effect to be of order 2.5 pixels over the 1000-pixel length 
of the CCD, and the arc lines are tilted by approximately 2 pixels. The distortion in each object 
and standard star frame was found by tracing the continuum along the dispersion direction. The 
distortion in the slit direction was similarly found by tracing the lines in the arc frames. Two 
dimensional fits to the wavelength and geometrical correction as a function of position was done 
with the IRAF task fitcoords. The two transformations were inverted to map the wavelength as a 
linear function along lines of the frame, and the position along the slit as a linear function along 
the columns. The subtraction of the background sky contribution could then be performed, having 
first normalised the spectra by fitting a 2nd order polynomial to the background in each column. 

All frames were divided by a smoothed spectrum of a featureless DC white dwarf star in order 
to remove the effects of the dichroic transmission variability, which because of interference effects is 
a complex oscillatory function of wavelength. The sensitivity function derived from the measured 
standard stars can then be fit by a low order polynomial. The spectrophotometric standards L745- 
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46A and W485A were used to calibrate the absolute fluxes of the spectra on the first night, while 
only L754-46A was observed on the second night. The nights appeared to be photometric, with 
atmospheric transmission varying less than 10%. No attempt was made to remove the atmospheric 
molecular absorption bands. 

An error frame was calculated for each 2-dimensional spectrum by taking the square root of 
the counts after the instrumental corrections and propagating the error through the reductions. 

Since the observations were made without the aid of an auto-guider, and the slit positions 
were set up by offsetting from the nuclear peak of light as it appeared in the TV display, the slit 
positions were not very accurately controlled. More accurate slit positions were found by matching 
the intensity profile of the brightest line in the spectrum [O III] A5007 to the intensity profile 
generated from an narrow band image of the galaxy in this same emission line. This technique 
involves firstly rotating and binning the image to the same pixel scale and spatial orientation as 
the 2-d spectra ((X'89 pixel" 1 , PA 30°). Simulated intensity profiles for slits positioned at Of! 3 steps 
across the ENLR were then generated from the image, and were compared to the observed slit 
profiles. The best match in each case was found by cross-correlating the observed intensity profile 
with the set of simulated slit profiles. The IRAF task fxcor was used in its pixel correlation mode for 
this purpose. The correlation function calculated for each simulated slit with the observed profile 
is a measure of how well the simulated profile matches the observed profile, and gives the relative 
shift in pixels between the two profiles. The height of the correlation functions provide three or four 
best match candidates, which are then inspected manually and a best match chosen. For all but 
slits -5 and 3 the best match was obvious upon inspection so in general we regard the slit positions 
as known to within ff.'3. The derived slit positions are shown overplotted on the [O III] A5007 image 
contours in Figure 2. In the remainder of the paper we refer all spatial positions in the ENLR of 
NGC 2992 to the coordinate system shown in this figure. 

The relative and absolute flux calibrations of all slits was assessed by extracting the flux of 
[O III] A5007 from the calibrated narrow band image at the x, y coordinates of the spectra and 
comparing with the line flux measured from the spectra. The spectra and image fluxes match 
within 20% at most x, y positions. 

One dimensional spectra were extracted from the long slit frames by summing every three 
pixels along the slit. The y coordinate of the extraction is known from the y position of the best 
match simulated slit, and the x coordinates are determined from the relative shift given by the 
cross-correlation. This results in a grid of spectra covering the ENLR of NGC 2992 with a spacing 
of approximately 2" in the y direction, and closer spacing in the x direction. To improve the signal- 
to-noise ratio in the faint, outer regions of the ENLR we have co-added up to 8 spectra into a more 
coarse grid. 

Average spectra of the NW cone, SE cone and inter-cone regions of the ENLR are shown in 
Figures 3a, 3b and 3c. A sky feature is seen at 6860 A and the red end of the spectra are not 
shown as that region is contaminated with sky features. 
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2.1. Emission line measurements 

The fluxes, wavelengths and widths of the bright emission lines were measured in each grid 
point spectrum by fitting Gaussians to the line profiles with a linear fit to the continuum in the 
wavelength region of the line. We find that approximating the continuum as a linear function does 
not significantly affect the measurement of the bright lines. The bright emission line measurements 
were done with the IRAF STSDAS task specfit. For the [S II], [N II] and [O III] doublets the 
FWHM and wavelength ratio of the lines were fixed in the fitting procedure, la uncertainties in 
the flux, wavelength and FWHM measurements were calculated from the error spectrum derived 
from the error frames described above. 



2.2. Subtraction of the stellar population and faint line measurements 

At all locations in the ENLR the emission line spectrum is superimposed on an underlying 
galaxy stellar spectrum. Figure 3b shows that the galaxy contribution is much stronger on the SE 
side of the galaxy. The main stellar features are readily identified as the Ca H and K bands, the 
G band of C-H at 4290-4313 A, H/3 absorption, Mg b at around 5170 A and Na I 5890-5896 A. 
We find that these stellar features and the general shape of the underlying galaxy spectrum can be 
matched by the spectrum of the E4 peculiar galaxy NGC 1316 ( Pickle^ 1985) with some correction 
for reddening. 

Subtraction of the stellar population is very important for the detection and measurement of 



faint emission lines. ptorchi-Bergmann et ajj (1996) emphasize the importance of subtracting stellar 



features to obtain accurate measurements of the H7, [O III] A4363, and He II A4686 lines in ENLR 



spectra. Using Storchi-Bergmann et al.| (1996) as a guide for matching the stellar features we have 



fitted the observed spectra of the NW and SE cones, using the NGC 1316 galaxy template to account 
for the underlying galaxy continuum. The fit was done over the wavelength range 4200-5300 A, 
with the IRAF specfit package. In the fitting procedure we constrain the redshift of the template 
to a range close to that determined for the emission lines in the spectrum. A reddening correction 
is applied in the fitting procedure using the Cardelli, Clayton, fc Mathis| (1989) reddening curve 



and allowing E(B-V) to be a free parameter. The emission lines H7, [O III] A4363, He II A4686, H/?, 
the [O III] AA4959,5007 doublet and [N I] A5200 were included in the fit as Gaussian components. 
The widths of the H7 and H/3 lines were linked to have the same value in the fit, as were the 
widths of the [O III] AA4959,5007 doublet lines. [O III] A4363 was not linked to the other [O III] 
lines because the asymmetry observed in the brighter [O III] lines is not seen in the the fainter 
[O III] A4363 line. An error spectrum derived from the original photon statistics of the CCD 
observations was used in the fitting procedure and allows the calculation of uncertainties in all the 
fit parameters. The la uncertainties in the NW and SE cone spectra are approximately 0.15 and 
0.4 xlO -16 ergs s _1 cm -2 A -1 respectively over the wavelength range of the fit. 



The NW cone and SE cone spectra are both well fit by this procedure, with the template 
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galaxy component matching the wavelengths and depths of the stellar absorption features. The 
underlying galaxy spectra in the NW cone required an E(B-V) of 0.16 and the SE cone required an 
E(B-V) of 0.17. Figures 4a and 4b show the result of the fits for the NW and SE cones respectively 
where we show the galaxy template component overplotted on the observed spectra. For each case 
we have also subtracted the fitted galaxy template spectrum to reveal the spectrum of the emission 
line gas. A small amount of residual structure is seen the spectrum after template subtraction, 
however this structure is within the la noise level. The template subtracted spectrum of the NW 
cone clearly shows the H7, [O III] A4363 and He II A4686 lines. In the SE cone H7 and He II A4686 
are clearly seen, but [O III] A4363 is only just above the noise level resulting in a relatively uncertain 
measurement of the [O III] A4363 flux for the SE cone. 



2.3. Reddening correction 

NGC 2992 has a prominent dust lane and the observed spectra show that the reddening varies 
as a function of position in the ENLR, and that the reddening derived from the emission lines is 
higher than for the underlying galaxy spectrum. The most severely reddened spectra occur in the 
regions associated with the dust lane, with the Ha/H/3 ratio observed as high as 10. For each grid 
point spectrum where both Ha and H/3 are detected we derive a value of the reddening from the 



Ha/H/3 ratio, using the reddening curve of Cardelli, Clayton, fe Mathis (1989) and assuming an 



intrinsic Ha/H/3 ratio of 3.1. The reddening within the dust lane thus corresponds to Ay=3.6 and 
decreases outwards along the y direction away from the dust lane. The average reddening values 
in the NW and SE cones are Ay=0.86 and 1.8 respectively. With the orientation of NGC 2992 
described in section 1, we might expect the NW cone, seen through the disk of the galaxy to have 
the greater reddening. The disturbed nature of the dust lane may account for this discrepancy. 

The emission line measurements for the average spectra of the NW and SE cone, and inter-cone 
regions are given in Table 1 along with their reddening corrected values. In each case the line fluxes 
are given relative to H/5=100. 



3. RESULTS 

The spectra of the ENLR of NGC 2992 are typical of Seyfert 2 type galaxies with strong lines 
of [O III] A5007, [O II] A3727, Ha and also the low ionization lines [S II] AA6717,6731, [N II] A6583 
and [O I] A6300. Highly ionized Iron lines (Fe VII, Fe IX, Fe X) sometimes seen in Seyfert galaxies 
are not detected in our spectra. [Ne V] A3426 is detected at the far blue end of the wavelength 
range, and the [S III] AA9069,9532 lines are detected but are strongly affected by night sky features. 
We also detect a blue asymmetry in the line profiles of the [O III] AA4959,5007 lines in the SE cone. 
The fainter lines of [O III] A4363, H7 and He II A4686 are also detected in co-added spectra of the 
NW and SE cones. 
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3.1. Spatial variation of the spectra 

A striking feature of this set of spectra is the strong variation of the line ratios as a function of 
position in the ENLR. Figure 5a shows how the [O III] A5007/H/3 and [N II] A6583/Ha line ratios 
separate the spectra of the ENLR into three groups that correspond to different regions within 
the ENLR. The open circles indicate spectra extracted from the NW cone of the ENLR, the filled 
squares are for points within the SE cone, and the open triangles correspond to spectra extracted 
from regions between the two cones, and from the outer edges of the two cones. 

This effect is shown even more clearly in Figure 5b where we have used the ratio of [O III] A5007 
to Ha instead of H/3. Using the brighter Ha line allows us to include spectra down to a fainter 
brightness limit, thus covering more of the ENLR, especially those fainter regions between the cones. 
The disadvantage of using the [O III] A5007/Ha ratio is that the reddening correction is greater than 
the correction to [O III] A5007/H/3. This larger reddening correction with its greater uncertainty 
causes more scatter of the points in Figure 5b. For the faintest regions where H/3 is not detected, we 
have applied a value for the reddening which has been extrapolated from adjacent points at which 
H/3 is measured. This procedure is most uncertain for the spectra extracted from points between 
the cones and within the dust lane of the galaxy. For these points we apply the average reddening 
of the dust lane points which is equivalent to assuming an Ha/H/3 ratio of 6 before reddening 
correction. Despite these uncertainties in the reddening correction for the [O III] A5007/Ha ratio, 
the inferred spatial variation of the line ratios remains similar to those found in Figure 5a. 

Figure 6 shows the spatial locations of the spectra used in Figures 5a and 5b. We shall 
maintain the distinction between these three spatial regions by using the same plotting point styles 
throughout the rest of this paper. 

3.2. [N II] /Ha variation 

Figures 5a and 5b show that the the NW and SE cones have significantly different [N II] A6583/Ha 
ratios. The [N II] A6583/Ha ratio of the average spectrum of the NW cone is 0.46 compared to 
0.76 for the SE cone. 

Since the SE cone has a stronger underlying galaxy continuum we have investigated the pos- 
sibility that the variation in [N II] A6583/Ha between the two cones is due to a greater amount 
of Ha absorption in the SE cone. To test this possibility we use the fit to the blue region of the 
underlying stellar distribution that was discussed in the previous section, and apply it to the spec- 
tral range around the Ha and [N II] lines. In the spectral region around Ha it is more difficult to 
make an independent fit for the underlying galaxy continuum. We therefore fit the galaxy template 
in the 6450—6700 A spectral range using the same multiplicative scale and reddening parameters 
derived for the blue region. The Ha line and each of the [N II] AA6548,6583 emission lines are 
fitted with single component Gaussian profiles with the fit constrained so that the [N II] doublet 
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line wavelengths were at the known ratio. For the fit the flux ratio [N II] A6583/A6548 was fixed 
at 3.0, and both [N II] components were constrained to have the same line widths. 

The results of the fits to the NW cone and SE cone and dust-lane region spectra over the 
6450—6700 A range are shown in Figure 7. In each panel we show the total fit overlaid on the data 
and also the result of subtracting the template from the data. The lighter weight curves in Figure 7 
show the fitted galaxy template and the individual Gaussian component fits to the Ha and [N II] 
emission lines. 

We obtain an excellent fit to the spectrum in the NW cone. Since the stellar contribution 
is faint in this cone there is only a small amount of Ha absorption which is barely discernable in 
the plot, and is insignificant relative to the emission line fluxes. The fit to the SE cone (middle 
panel of Figure 7) demonstrates that stellar Ha absorption is more important on the SE side of 
the ENLR. This has the effect of increasing the Ha flux measurement by < 5% and decreases the 
[N II] A6583/Ha ratio by < 7% compared to the values obtained by using a linear continuum level 
in the measurement of these lines. This is not large enough to make a significant difference to the 
[N II] A6583/Ha ratio. We have thus established that the variation of the ratio between the NW 
and SE cones cannot be attributed to Ha absorption, and must therefore be due to real physical 
differences between the excitation conditions in the cones. 

The slight excess shown by the fit at the wavelength of the [N II] A6548 line is most likely 
due to a blue asymmetry in the Ha line similar to the asymmetry in the [O III] A5007 line seen 
at the same location. While it is possible to add another component of Ha to obtain a better fit 
to the spectrum, we would also need to include similarly blue shifted [N II] components, and this 
data simply do not have the resolution or signal to noise required to properly constrain the relative 
fluxes or line widths of such additional components. 

The lower panel of Figure 7 shows the fit to the brightest point in the dust-lane region of the 
ENLR. Here Ha and the [N II] lines are well fitted by single component Gaussians, and yield a 
[N II] A6583/Ha ratio of 0.59 ±0.01 which is intermediate to the values measured in the NW and 
SE cones. 



4. Physical properties of the emission line gas 
4.1. Temperature 

The [O III] A4363/[0 III] A5007 line ratio is the well known temperature diagnostic often re- 
ferred to as Rom- We estimate the temperature from Rom using the IRAF STSDAS nebular 
package which draws on the atomic data of Mendozaj (1983) for the +2 ion. An assumed density of 



100 cm" 3 was used in the temperature calculation, and we note that the calculation is relatively in- 
sensitive to density. The average spectrum of the NW cone has Ro///=0.019 ±0.004 corresponding 
to T e =14780lggQ°. The average spectrum of the SE cone has Rom=0.025 ±0.009 corresponding 
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to T e =16990_23ooi the greater uncertainties in the SE cone are due to the large uncertainty in the 
[O III] A4363 measurement. These temperatures are lower than the value of 19900 K obtained by 
purret and Bergeron| (1988) for the nuclear region. It is well known that temperature fluctuations in 
the ionized region tend to drive the value of Rom up, so the measured temperatures are probably 
somewhat higher than the mean temperature characterising the +2 zone. 



4.2. Density 

The ratio of [S II] A6717/[S II] A6731 is a density diagnostic sensitive in the range 100—10000 cm -3 . 
Density estimates from the [S II] ratios were calculated with the IRAF STSDAS nebular package 
which draws on the atomic data of |Cai Pradhan (1993) for the S +1 ion. Figure 8 shows these 



densities plotted against the projected radius of the spatial locations within the ENLR. The average 
error bar on the density measurements is approximately 80 cm -3 . We find that the [S II] ratios are 
close to the lower limit of the diagnostic range, and that some points have [S II] ratios which can 
only be used as an upper limit for the density. The central regions have the higher density values. 
Although Figure 8 shows a general trend of decreasing density with distance from the nucleus, the 
densities in the cones (when these can be measured) remain relatively constant with distance. 



4.3. Ionization Parameter 

I 1 

Penston et al. (1990) demonstrated that the [O IIl]/[0 II] line ratio is proportional to the 
ionization parameter U for radiation-bounded photoionized clouds irrespective of the shape of the 
ionizing continuum. Following Penston et al. (1990) we estimate the ionization parameter by the 
relationship: 



U = 1.82 x 10~ 3 [O III] A5007 / [O II] A3727 (1) 

where we use the reddening corrected ratio of [O III] A5007 / [O II] A3727. Figure 9 shows 
the values of U plotted against the projected radius. Interpreted literally, this shows that the NW 
cone has slightly higher mean U than the SE cone, while the dust lane and outer regions of the 
cone have the lowest ionization parameter. 



5. LINE RATIO DIAGRAMS 

The utility of line diagnostic diagrams consisting of selected pairs of emission line ratios is well 
established. With such diagrams, different excitation mechanisms can be distinguished, chemical 
abundances may be estimated (in special circumstances), and physical parameters may be derived. 
The most commonly used diagrams involve ratios of bright, easily measured lines ( [Baldwin, Phillips] 
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&; Terlevich 1981), and the best of these restrict the choice of line ratios to emission lines that are 
close together in wavelength, avoiding the necessity to determine reddening corrections with great 
accuracy ( Veilleux and Osterbrock| 1987). 

Our main objective in this section is to analyse the excitation mechanisms responsible for the 
emission lines. To do this we compare the observed line ratios with grids of shock and photoion- 
ization model predictions. 

The models used here are: 



The MAPPINGS II shock and shock+precursor models of Dopita Sutherland! (1996). These 



models are parameterised by the shock velocity V^ oc fc=200— 500 km s _1 , and the magnetic 
parameter B/n 1 ^ 2 = 0—4 /xG cm 3//2 . In these models, EUV photons produced within the 
cooling region of the shock, diffuse both upstream and downstream. Those that go down- 
stream are trapped by the shocked gas and produce a low-ionization parameter photoionized 
region near the point where hydrogen recombines. The combined spectrum of the cooling 
and photoionized gas we refer to as the "shock only" model. The (undoubtedly important) 
effect of thermal instabilities are not accounted for in this model. The photons which pass 
upstream produce a photoionized region of high ionization parameter and excitation in the 
pre-shock medium. The importance of this region depends upon whether this medium is able 
to absorb the ionizing photons. Shock models which include the contribution of both the 
shocked gas and its ionized precursor are referred to as "shock+precursor" . 



The power-law photoionization models presented in Allen, Dopita fc Tsvctanovj (1998). These 



are simplified photoionization models calculated with a ionizing spectrum of the form oc v a , 
with a=-l, -1.4 and -2 and a range of ionization parameters (10~ 4 < U < 1), and densities 
(100—1000 cm -3 ). These models were calculated with the same MAPPINGS II code as the 
shock models described above. 



The A M /i sequence of Binette, Wilson, & Storchi-Bergman (1996). This is a sequence of 



photoionization models which include both matter bounded and ionization bounded clouds, 
such that the ionizing spectrum incident on the ionization bounded clouds is filtered by the 



matter bounded component. Binette, Wilson, &: Storchi-Bergman| (1996) find that this two 



component model can account for some of the discrepancies of single component photoion- 
ization models. This sequence is parametrised by A^/j which is the solid angle ratio of the 
matter bounded clouds to the ionization bounded clouds. 



In plotting the observed points on the resulting line ratio diagrams, we maintain the same 
plotting point styles as used in section 3 to distinguish between the NW cone (open circles), SE 
cone (closed squares) and inter-cone (open triangles) regions. In each diagram we show the average 
error bar of the observed line ratio, and also provide an arrow indicating the direction of reddening 
correction. The shock model grids are labelled with the shock velocity and magnetic parameter. 
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The power-law photoionization U sequences are plotted for n = 100 cm -3 and n = 1000 cm -3 with 
grid lines of constant U every 0.25 dex. The A M /j sequence is labelled with A M /j and tick marks 
at intervals of 0.2 dex. 



5.1. The Bright-line Diagnostic Plots 

In Figures 10 through 16 we show the bright-line diagnostic plots which are generally used 
for determining excitation conditions in AGN. On these diagrams, we have sufficient observational 
points to analyse the detailed spatial structure of the ionization cones and the inter-cone region. 
In all of these diagrams the line ratios fall within limits that characterise Seyfert galaxies, in fact 
the points from this single object cover almost completely the regions of this diagram filled by 
samples of Seyfert galaxies ( |Veilleux and Osterbrock" 1987), indicating that the NLR of NGC 2992 



is characterized by a wide range of physical parameters. This is also reflected in the way that the 
points fall over a large range of the parameter space covered by the various models. 

In section 3 we have already shown how the [O III] A5007/H/3, and [O III] A5007/Ha versus 
[N II] A6583/Ha diagrams separate the NLR of NGC 2992 into 3 distinct spatial regions. These 
diagrams are repeated here (Figures 10 and 11) with the model grids overlaid. Similar behaviour 
is clearly shown in other excitation-dependent diagnostic plots such as [O III] A5007/H/3 versus 
[S II] AA6717,6731/Ha (Figure 12), or [O III] A5007/H/3 versus [O II] A3727/[0 III] A5007 (Fig- 
ure 13), but is not apparent when the ratios of low-excitation lines with hydrogen are compared, 
such as the [N II] A6583/Ha and [S II] AA6717,6731/Ha line ratios (Figure 14). In these diagrams 
the points from the different regions tend to cluster more tightly, and are not clearly separated as 
they are in Figures 10 and 11. 

We will now consider how well the various theoretical sequences match (or do not match) the 
observations. First, the shock-only grid fails to match any of the observational data, for any of 
these diagnostic plots. Thus, if the ENLR gas is shock-excited, the shocks must be occuring in a 
matter-rich environment where the EUV photons produced by the shocks can be absorbed. We 
will not discuss the shock-only grids further. 

The shock -l-precursor grids do much better. In all of the diagrams, with the exception of those 
which involve ratios of the [O I] lines, such as [O I] A6300/[O III] A5007 (Figure 15), the observations 
define shock velocities which are reasonably consistent from one diagnostic diagram to another. For 
example, in Figures 10a, 11a and 12a the points from the NW and SE cones fall within the high 
shock velocity region (350—500 km s _1 ) of the shock -l-precursor model grid. The inter-cone points 
show a large range in [O III] A5007/H/3 which extends these points along lines of decreasing shock 
velocity, down to as low as 200 km s _1 with approximately constant magnetic parameter in the 
shock -l-precursor model grid. However, the different diagnostic plots do not agree very well as to 
the value of the magnetic parameter, other than broadly excluding the possibility that this is less 
than about B/n 1 / 2 ~ 1/zG cm 3 / 2 . However, such scatter as is observed in this parameter between 
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the different diagrams can probably be ascribed to differences between the abundance set used in 
the models, and those that are applicable to NGC 2992, as well as to collisional de-excitation effects 
in the [O II] and [S II] lines, which are not modelled in these (low density limit) shock models. 

Turning now to the simple photoionization models. In each of the diagnostic plots, the pho- 
toionization models overlap the same region of parameter space in this diagnostic diagram as 
the shock+precursor models, and so the observed points also fall within the photoionization U- 
sequences. However, different regions tend to be characterized by different ionization parameters 
and spectral index. For example, on Figure 10b, the NW cone points are within the a = —1.4 
models with ionization parameter 1CT 2 ' 5 < U < 10~ 2 ' 25 . The points from the SE cone show more 
scatter, but indicate an average ionization parameter which is lower than in the NW cone. Most of 
the SE cone points fall between the a = —1.4 and a = —1.0 sequences which may indicate a higher 
spectral index in the NW cone than in the SE cone. There is however a degree of degeneracy in the 
photoionization models such that increase in density or an increase in spectral index have a similar 
effect on the predicted line ratios. 

The spread of the inter-cone points may therefore be partially explained by a wider range in 
density than we computed in the photoionization models. To investigate this possibility we extended 
the photoionization grid by calculating models with densities down to 10~ 2 cm -3 , for ionization 
parameters U = 10 -2 and 10~ 2 ' 75 , and spectral index a = —1.4. The low density limit occurs for 
density n e ~1.0 cm~ 3 , below which [O III] A5007/H/3 remains constant at ~ 10 06 for U = 10~ 2 
and ~ 10 09 for U = 10 -2 75 . However, the observed points extend over a wider range down to 
[O III] A5007/H/3 ~ lO 3 , making it difficult to account for the observations with a photoionization 
models with a single value of the spectral index a such as would be required in a self-consistent 
photoionization model of the whole NLR. 

Different diagnostic plots yield different estimates for both a and for U. For example, the NW 
cone points that fall within the a = —1.4 [/-sequence in Figures 10b and lib, are shifted closer 
to the a = — 1 [/-sequences in Figure 12b, with a lower value of U. The shift of the points with 
respect to the photoionization [/-sequences is significant, being approximately 4 times the average 
la error bar of the [S II] AA6717,6731/Ha observations. For the [S III] AA9069,9532/Ha versus 
[S II] AA6717,6731/Ha plot (Figure 16), a ~ —2 is indicated. In general, most of the plots arc 
consistent with U ~ 10 -2 - 5 for the bright parts of the cones, but those involving [N II] A6583 /Ha 
versus [S II] AA6717,6731/Ha, and [S III] AA9069,9532/Ha versus [S II] AA6717,6731/Ha are more 
consistent with U ~ 10~ 4 . 

We can therefore conclude that not only are the observations inconsistent with a single value of 
a within individual line ratio diagrams, but also that different diagnostic diagrams do not furnish 
consistent values of a between one another. Since photoionization models should be characterised 
by a unique value of a determined by the intrinsic nuclear spectrum (although U may change from 
one region to another), these observations provide a serious challenge to photoionization modellers. 

By contrast, the shock+precursor models indicate values of the shock velocity which are con- 
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sistent from one diagram to another; of order 350—500 km s _1 in the cones, falling to as low as 
200 km s _1 in the inter-cone region. In the shock+precursor picture, the apparent changes in both a 
and U according to the photoionization interpretation is the result of changes in the spectral energy 
distribution and the total energy flux in the ionizing radiation as a function of shock velocity. 

These problems with photoionization models are not resolved by use of the A M /j sequence. 
Indeed, the fact that A M /j sequences are very similar to [/-sequences, when plotted on the tra- 
ditional line ratio diagrams was shown in Binette, Wilson, fc Storchi-Bergman| (1996). This is 



apparent in Figures 10b, and lib where the A M /j sequence is very close to the higher density 
a=-lA [/-sequence. With respect to the Am/j sequence, the observed points from the SE and 
NW cones shown in Figures 10b and lib cluster around the sequence for values of A^/i in the 
range 1<^ A M /j <^6. The SE cone points tend to have [O III] A5007/H/3 higher than the A m jj 
sequence prediction, while the NW cone points have [O III] A5007/H/3 slightly lower than the 
Am/ / sequence, and favour higher values of A M /j. As with the [/-sequences, the inter-cone points 
tend to diverge from the Aj^n sequence. In general the A^/i sequences can neither explain the 
range or the direction of scatter of the observations in NGC 2992. In addition, in some cases such 
as in the [N II] A6583/Ha versus [S II] AA6717,6731/Ha and the [S III] AA9069,9532/Ha versus 
[S II] AA6717,6731/Ha diagrams (Figures 14 and 16 respectively), the A M /j sequence is offset from 
the observations by an amount which is larger than could be explained by abundance variations. 

The result of comparing the ratios of the bright lines mapped over the NLR of NGC 2992 can 
be summarized as follows: The ratios of [O III], [N II], [S II] and [S III] generally fall within the 
parameter space covered by both the shock+precursor, and photoionization [/-sequence models. 
The ratios of these lines fall consistently in the same location of the diagnostic diagrams relative 
to the shock+precursor grid, but do not provide consistency with respect to the photoionization 
models. [O II] appears to be slightly overestimated in both shock+precursor and photoionization 
models, and [O I] is predicted to be too strong in all the models. 



5.2. Faint line diagnostic diagrams 

The lines of [Ne V] A3426, He II A4686 and [O III] A4363 also are important diagnostics, but are 
relatively faint lines so we are not able to spatially map them over the entire ENLR of NGC 2992. 
In the following three diagrams we compare the model predictions with the observations of these 
lines in the average spectra of the NW and SE cones. 

Figure 17 shows [S II] AA6717,6731/Ha versus He II A4686/H/3. In this diagram the shock- 
only and shock+precursor models overlap, so we show the shock-only model grid with dashed 
lines. It is clear from this diagram that the He II A4686/H/3 ratio is underestimated in the 
shock+precursor models but the photoionization [/-sequence shown in Figure 17b can account 
for the large He II A4686/H/3 ratio with a = — 1 and ionization parameter 10 -3 < U < 10~ 2 . 



The line ratios used in Figure 17 are identified by Binette, Wilson, fc Storchi-Bergmanj (1996) 



-17- 



as one of the combinations for which A^/i and [/-sequences show markedly different behaviour. 
This is apparent in Figure 17b where the A M /j sequence cuts across the [/-sequences curves. The 
data points for the NW and SE cones fall relatively close to the A m jj sequence for Aj^n 

Figure 18 shows He II A4686/H/3 with [O II] A3727/[0 III] A5007, where we have plotted 
the shock models, photoionization models and A^/i sequence on the one diagram. The relative 
positions of the various grids in this diagram are similar to those Figure 17, and we find that the 
points fall in similar positions with respect to the A m jj sequence and [/-sequences, but fall well off 
the shock and shock -l-precursor grids. 

Figure 19 combines the ratio of the high excitation [Ne V] A3426 line to [Ne III] A3869 versus 
the [O III] A5007/H/3 ratio. The data points from the SE and NW cones fall approximately 3<r to 
the right of the shock -l-precursor models, but do fall within the a = — 1 photoionization model. The 
inferred ionization parameter from the position of the data points on the [/-sequences is much higher 
than inferred from the other diagrams, however the fact that the U-sequences fold over themselves 
in this region of the diagram (causing some regions to have double valued solutions within the 
photoionization models) makes this inference very uncertain. The Aj^n sequence again shows 
different behaviour to the [/-sequences and the data points fall near the Am/j ~0.6 prediction. 

Shock models of this simple class are known not to provide a good description of the very 
high-excitation line spectrum. This includes the He II and [Ne V] lines as well as higher excitation 
species such as [Fe VII]. These lines are better predicted by models which include a component 
of optically thin gas with high ionization parameter. The Aj^n sequence is an example this. 
The two components used in the A M /j models can be considered to be a simple approximation 
of an inhomogeneous medium which may contain regions of different optical thickness to ionizing 
radiation. The success of the A^/i sequence in predicting the line ratios as shown in the faint- 
line diagnostic diagrams may attest to the presence of inhomogeneities in the ionized clouds. The 
ionized precursor gas in our shock -l-precursor models is homogeneous and is thus not expected to 
reproduce these lines accurately. Nevertheless, Dopitaj (1997) has shown that an optically thin, 



high ionization parameter component which emits high ionization lines, can be produced in the 
context of more complex shock models associated with shocked entrainment layers at the boundary 
of ionization cones. 

It may also be possible for combined shock and central source photoionization models to 
provide a better match to the observed line ratios. Such models are however difficult to constrain 
and we find no obvious spatial or spectral separation of shock and central source components in our 
data, so we have not attempted a combined model. Nevertheless it is clear that these diagnostics 
involving faint high excitation lines will be very useful in developing more sophisticated models of 
shocks, photoionization or combinations of shocks and central source phtoionization. 
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5.3. The temperature sensitive diagnostic diagram 

A long-standing problem in AGN modelling has been the so-called temperature problem. 
Briefly the nature of this problem is that simple photoionization models give low temperatures 
(less than 15000 K), which disagrees with observations of [O III] A4363/[0 III] A5007 in NLRs of 
some active galaxies, for which the inferred electron temperatures range up to ~22000 K ([Tad 



hunter, Robinson fc Morganti| 1989). Shocks can better account for the high temperatures because 



the [O III] A4363 emission from shocks is dominated by the 40000-20000 K region of the cooling zone 



behind the shock front (Dopita & Sutherland 1995). Recently Binette, Wilson, & Storchi-Bergman 
(1996) have shown the A M /j models can also predict high values of [O III] A4363/[0 III] A5007. Also 
some of the earlier [O III] A4363 line measurements may be overestimated due to poor subtraction 
of the underlying galaxy spectrum. We now compare the temperature sensitive line ratios observed 
in NGC 2992 to the model predictions. 

Figure 20 employs the temperature sensitive ratio of [O III] A4363/[0 III] A5007 on the y- 
axis, versus He II A4686/H/3 on the x-axis. Due to the different temperature regimes of the shock 
and photoionization models, the grids are widely spaced in this diagram. A temperature axis 
is plotted on the left side of the diagram showing the inferred temperatures calculated from the 
[O III] A4363/[0 III] A5007 ratio in the same manner as discussed in section 4.1. 

The large [O III] A4363/[0 III] A5007 values predicted by the shock-only model reflect the 
high temperatures that occur in the cooling region behind the shock front. The photoionization 
model grids display the well known result that [/-sequences generally have [O III] temperatures less 
than 15000 K. The highest [O III] temperatures obtained in photoionization models occur for the 
a = — 1 model with ionization parameter U = 0.1, with the A^/i able to produce a similarly high 
[O III] temperature. As shown in section 4.1 the average spectra of the NW and SE cones have 
[O III] temperatures of 14780 K and 16990 K respectively. The location of the data points on this 
temperature axis are consistent within the error bars with the [O III] temperatures predicted by 
the shock models, the A M /j sequence and also the a=-l, [7=0.1 photoionization models. 



Binette, Wilson, & Storchi-Bergman (1996) show that high temperatures are obtained in the 
Am 1 1 sequence because the MB component is chosen to be sufficiently thin that the temperature 
is governed by the photoionization of He + , and that the +2 zone is truncated before the internal 
temperature of that zone falls enough to decrease the [O III] A4363/[0 III] A5007 ratio. The A m jj 
sequences described by Binette, Wilson, fc Storchi-Bergman| (1996) (and displayed in Figure 12 



of [Binette, Wilson, k Storchi-Bergmanl 1996) can reach [O III] A4363/[0 III] A5007 ~0.019 for 



an incident a = —1.3 ionizing spectrum (as plotted in Figure 20), up to ~0.025 for an incident 
ionizing spectrum with a = —1.1, and as high as ~0.048 when the metallicity of the gas is decreased 
to Z=0.2. As such these A M /j sequences are able to reproduce the relatively high temperatures 
observed in the NW and SE cones of NGC 2992. 
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6. CENTRAL SOURCE TEST 
6.1. Determination of Q from [S II] density. 

I I 

Metz et al. (1996) describe a method to test whether the line emitting gas in the NLRs of 
Seyfert galaxies is photoionized by a central source. The test uses the definition of the ionization 
parameter for radiation from a central source incident on an emission line cloud: 



U = ~V > ( T m ( 2 ) 
n e r z c n e {r/ cos c/J 

where q is the number of ionizing photons per unit time per steradian, n e is the electron density, 
and r is the distance of the emission line cloud from the source. 6 is the angle between a line from 
the source to the cloud and the plane of the sky. Using [O III] A5007/[O II] A3727 to estimate U, 
[S II] A6717/A6731 to measure the density n e and estimating r, allows q to be determined. If the 
assumption of a central source is correct, then q should be independent of position in the NLR. 



Using long slit spectroscopy Metz et ah] (1996) applied this technique to the NLR of three nearby 



Seyfert galaxies including NGC 2992. In the case of NGC 2992 their slit was aligned along the 
radio structure (PA = 26°). They find that q rises with increasing distance from the nucleus in all 
three NLR. 

Here we perform a similar test on NGC 2992 using our data which has more complete spatial 
coverage of the NLR than the single slit in Metz et al.| (1996). We confirm the result that q 



measured in this way does rise with increasing distance from the nucleus. However we find that 
the [S II] A6717/A6731 density diagnostic is not sensitive in the density range required to measure 
q in the NLR of NGC 2992, and thus demonstrate that the test is invalid. 

To perform the "g-test" with our spectra of NGC 2992 we have considered only the regions 
within the bi-cone of opening angle 116° centred on the nucleus in which have well determined [S II] 
density measurements. U is determined by using the reddening corrected [O III] A5007/[O II] A3727 
ratio as described in section 4.3. The values of n e and U for points within the cone are shown in 
Figures 21a and 21b respectively. To estimate the de-projected distances of points in the bi-cone 
to the nucleus, we have considered all the points to lie within the geometric volume of the bi-cone. 
Narrow band images show that the bi-cone axis is aligned at roughly right angles to the plane of 
the galaxy. Here we assume that the bi-cone has the same inclination on the sky as the galaxy, 
and that the SE cone is the near one. Using this geometry we can then calculate two values for 
the de-projected radius of each point 0, which correspond to the near and far surfaces of the cones. 



2/(cos(V ) + i) + sin ip sin i) ± y 2 sin 2 ip sin 2 i — x 2 sin 2 i cos 2 (ip + i) — 2a; 2 cos( , + i) smip sin 3 i 

Tdeprojected = T7~i \ ^\ \ ~n 7~i \ -N '• I '• '• V J 

cos (ip + i) + 2 cos(V + i) sin tp sin i 
where i is the angle subtended by the cone axis to the plane of the sky, and tp is the half angle of the cone. 
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These two values represent the extrema for the de-projected radii of points within the cones, and 
we calculate q using both values, and likewise with the projected radius. 

The diagrams on the left hand side of Figure 22 show the values of Q = 4nq derived for; (i) 
de-projection onto the near side for both cones (top panel), (ii) de-projection onto the far side of the 
cones (middle panel) and (hi) no de-projection (bottom panel). In all cases Q is found to rise with 
increasing distance from the nucleus, and the result is relatively insensitive to the de-projection 
used. 

A rising Q is inconsistent with photoionization from the central nucleus, because in this case Q 



should have the same value wherever it is measured. Metz et al. (1996) rule out the possibility that 



this result is due to a projection effect because it would require 6 to increase towards the nucleus, 



which is inconsistent with the observed sharp edged shape of the cone. Metz et al. (1996) also 
note that a two-component ionizing continuum due to a mixture of ionization bounded and matter 
bounded clouds cannot account for the observed increase in Q. Although photoionizing shocks 



( Dopita fc Sutherland! 1995) provide a local source of UV which might explain the rising trend of 



Q, it is more likely that the apparent increase in Q is spurious, and is the result of overestimating 
the average density using the [S II] A6717/A6731 ratio. 

The [S II] A6717/ A6731 ratio is only sensitive as a density diagnostic over the range ~100— 10000 
cm -3 . The ratio tends to a limit of 1.4 for lower densities, and 0.44 for higher densities. In NGC 2992 
we find that the [S II] density for some locations near the nucleus is in the range 100—400 cm -3 , 
however for many locations the [S II] A6717/A6731 ratio lies near or below the lower limit of the 
usable range of this density diagnostic. If we take the value of Q as measured near the nucleus 
where the density is relatively well determined, we find Q = 10 52 ' 9 , in units of photons s _1 . This 
value is probably reliable, and reflects the true value for the central ionizing source. Using this value 
of Q, and the measurements of the ionization parameter and de-projected distances, we rearrange 
equation 2 to calculate the density required at each location, n e) g, to keep Q constant: 



log(n e ,Q) = log(52.9) - log(4vrc) - log U - 2 log r (4) 

The diagrams on the right hand side of Figure 22 show the results of this calculation, with 
log n et Q plotted against r for each of the de-projections. The dashed line indicates the lower useful 
limit for the [S II] density diagnostic. Most of the points indicate that n e ^Q is below this limit. 
Thus if Q is actually constant at 10 52 ' 9 photons s , we would not be able to measure it at large r 
because the density diagnostic is not sensitive in the range required. The g-test is therefore invalid 
in this case. However, if the local density is really as high as 200—300 cm" 3 in the regions at large r 
where we can measure it, then we would still require shocks if only to produce a local compression. 
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6.2. Determination of Q from the Ha luminosity 

In this section we compute the rate of hydrogen recombinations Q rec from the Ha luminosity 
and compare this to the rate of ionizing photons Q, derived from the emission line ratios above. 
In photoionization equilibrium the rate of hydrogen recombinations must equal the rate of ionizing 
photons. The rate of hydrogen recombinations Q rec can be derived from the Ha luminosity by 
comparing the effective recombination coefficient to produce Ha photons with the total (Case B) 
recombination coefficient for hydrogen: 

Q rec ~ 1.39 x 10 52 L 40 (Ha) photons s" 1 (5) 

where L/4o(Ha) is the Ha luminosity (in units of 10 40 ergs s -1 ). In the case where a central source 
is illuminating the NLR: 

4:7V 

Q rec ~ 1.39 x 10 52 L 40 (Ha) — photons s _i (6) 

ecu 

where e is the fraction of ionizing photons that are absorbed, and uj is the solid angle subtended 
by the Ha emitting region as seen from the nucleus. 

The reddening corrected Ha luminosities of the NW and SE cones are estimated from our 
Ha+[N II] narrow band image, corrected for reddening with a reddening map derived from the grid 
of Ha/H/5 ratios measured from the spectra. The reddening map was constructed by re-sampling 
the semi-regularly spaced grid of reddening constants from the spectra onto the image pixel grid 
with the IDL TRIGRID routine. The Ha+[N II] flux is summed over the specified region in the 
reddening corrected image, and a correction applied for the [N II] contribution using the average 
[N II] /Ha ratio of the region as measured from the spectra. This yields the total Ha flux which 
is then converted to a luminosity using the distance to NGC 2992 (30.8 Mpc). For each of the 
SE and NW cones, we have used only the Ha luminosity of the cones (corresponding to the same 
region over which Q is calculated in the previous section) and have specified u= 2.26 steradian, 
the solid angle of the cone subtended at the nucleus. The calculation of Q rec for the entire NLR 
region ignores the cone geometry and uses the total Ha luminosity and u = An. The dominant 
source of uncertainty in these calculations are errors in the reddening correction, and the implicit 
assumption that e is unity. This means that the Q rec values derived represent a lower limit estimate 
of the number of ionizing photons Q. 

The values of L4o(Ha) and Q rec computed for the NW and SE cones, and the NLR as a whole 
are shown in table 2^. Comparing Q rec derived from the Ha luminosity, and the Q derived from the 
forbidden emission line ratios we find that Qrec hes within the computed range of the Q values for 
all the regions considered. If the whole region was photoionized by a single central source then we 



4 Table 2 also contains L40QO III]) calculated from the reddening corrected [O III] narrow band image 
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should expect the values of Q rec to correspond to the value of Q as determined for the inner parts 
of the NLR in the previous section. However, our lower limits for Q rec are however consistently 
higher than the value of Q derived at the inner edge of the NLR (10 52 9 photons s _1 ). 

In conclusion, although the g-test in its original form is invalid, there does appear to be real 
evidence for a local source of ionizing photons with a hard spectrum within the cone. This can only 
come from photoionizing shocks. We note that a method along the lines of the g-test as described 
here, but devoid of the density diagnostic limitations, to quantify the relative contributions to the 
ionizing field from local and central sources. This would aid the development of combined shock 
and central source photoionization models. 

7. SHOCK ENERGETICS 

The emission line properties of the ENLR of NGC 2992 appear to be understandable in the 
context of shock excitation of the gas. In this section we compare the excitation of the gas to the 
kinematics of the outflow/wind, and also compare the energetics of the radio and X-ray plasma to 
the emission line gas. 

7.1. Excitation and Kinematics 

As described in section 1, the kinematics of the emission line gas shows a outflow velocity field 
superimposed on the underlying galaxian rotation. Figure 23a shows a map of the the velocity 
difference of the wind component to the galaxy rotation from our high resolution spectra of the 
[O III] line. The wind component is systematically blueshifted in the SE cone, and redshifted in 
the NW cone with velocities of up to ~200 km s _1 . A map of the F WHM of the wind component 
shown in Figure 23b reveals that the SE cone has the broader FWHM, and we consider this to be 
evidence of greater dynamical activity in the SE cone. 

If shocks are important for the excitation of the gas, we may expect the kinematical properties 
of the gas to in some way correlate with excitation. To compare the excitation of the gas with the 
kinematics we re-binned the high and low resolution spectra onto a common spatial grid covering 
the ENLR, and extracted the emission line ratios, and kinematical properties at each grid point. 
We now compare the [O III] A5007/H/3, and [N II] A6583/Ha ratios to the kinematics. 

In Figure 24a we plot the [O III] A5007/H/3 ratio, which is a measure of the excitation of 
the gas, versus the FWHM of the kinematical components. The kinematical component identified 
with the galaxy rotation are plotted as the solid circles, the main outflow components (as shown 
in Figure 23) are plotted as open squares and where the extra component is detected (at a few 
locations), this is plotted as open triangles. The uncertainty in the FWHM is generally less than 
100 km s _1 . 
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The outflow component points show a wide range in FWHM from 100 to 670 km s -1 , whereas 
the FWHM of the rotational component is much more quiescent, and is generally less than ~200 
km s _1 . The excitation of the rotational component covers the full range in [O III] A5007/H/3, while 
the wind components tend to show higher [O III] A5007/H/? ratios in locations where the lines are 
wider. No regions with line FWHM>250 km s _1 display low excitation. 

We also compare the data points on the excitation versus line width diagram to the theoretical 
curves predicted by the shock+precursor models. The solid curve assumes the FWHM is directly 
related to the shock velocity, and the dashed curve assumes that the FWHM is associated with 
the projection of randomly oriented shocks so that V s hock 

~FWHM/\/3. While the spread of 
the observed points is very suggestive of shock activity, it is not what would be expected for a 
simple shock moving into an undisturbed pre-shock medium. In this case, the regions of very high 
[O III] A5007/H/3 ratio would arise in the precursor, which would be kinematically undisturbed, 
and therefore there would be an inverse correlation of line width and excitation. However, if the 
shocks arise in a mass-entraining flow, the strong shocks are located near the outer boundary of 
the entrainment flow, and the photoionized precursor medium is the fast-moving material inside 
this ( Dopitaj 1997). This geometry would provide the direct correlation between line width and 



excitation observed. 

In Figure 24b we plot the [N II] A6583/Ha ratio versus the FWHM of the outflow components. 
As in Figure 24a the square points indicate the main outflow component, and the triangles indicate 
the extra component. Here we make the distinction between the cones by using filled symbols for 
the SE cone and open symbols for the NW cone. The cones are clearly separated in Figure 24b 
because the NW cone has lower [N II] A6583/Ha and lower FWHM, while the more dynamically 
active SE cone has higher [N II] A6583/Ha. The shock+precursor curves plotted on Figure 24b 
show that the models predict higher [N II] A6583/Ha for faster shock velocities. This suggests that 
the higher [N II] A6583/Ha in the SE cone is directly related to the greater dynamical activity in 
that cone. 

The difference in shock velocity between the cones inferred from these curves is however greater 
than that indicated by the [O III] A5007/H/5 ratio which indicates only slightly faster average shock 
velocity in the SE cone (see Figure 10a). As such, the increase in [N II] A6583/Ha is likely due to 
a combination of faster shock velocity and increased magnetic parameter, both of which may be 
expected in the more turbulent conditions within the SE cone. 



7.2. Thermal pressure energy source for the ENLR 

Another aspect of shock energetics is related to the question of whether there is sufficient 
energy in the radio plasma to power the optical emission. For Seyferts, classical equipartition 
arguments suggest that the pressure in the relativistic plasma is inadequate to drive high velocity 
shocks into the interstellar medium and produce the observed emission line flux (Wilson, Ward, & 
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Haniff 1988). Recently Bicknell et al.| (1998) have argued that the radio plasma may only trace a 



fraction of the available energy which exists in an entrained thermal phase, which is not detected 
in synchrotron emission, and that the NLR may well be powered by the energy and momentum 
flux originally carried by the radio-emitting jets. 



Briefly summarising the physics detailed in Bicknell et al. (1998): The total energy of a 
jet-fed lobe is f a d,FEt where Fe is the jet energy flux, t is the age of the lobe and f a( i ~ 0.5 
is a factor which accounts for energy lost in adiabatic expansion. Expansion of the lobe gives 
rise to narrow line emission via shocks and the jet energy flux Fe can be estimated from the 
[O III] A5007 luminosity, L([0 III])~ 0.01-Fg. The relativistic electrons in the lobe emit radio 
wavelength synchrotron radiation, and the ratio of the monochromatic radio power of the lobe (in 
synchrotron radiation) to the jet energy flux can be expressed as: 

k v ~ C(a)B^y"u-^^f e f ad t (7) 

where a is the power-law exponent of the electron energy distribution, C(a) depends on syn- 
chrotron parameters, and B is the randomly oriented magnetic field. f e is the fraction of the total 
energy in relativistic electrons and positrons, and is thus a measure of non-thermal to thermal 
composition of the plasma. We now use the relationships detailed in Bicknell et alj (1998) to 



investigate the possibility of a thermal energy source for the ENLR of NGC 2992. 

Our reddening corrected [O III] A5007 luminosity gives an estimate of Fe = 9.6 x 10 ergs s , 



and the 335 mJy flux density detected in the Fleurs map of |Ward et ag (1980) gives a radio power 



at 1.4 GHz of Pi. 4 =3.8 xl0 2y ergs s" 1 Hz" 1 . This yields k 1a = 4 x 10 -15 Hz" 1 , which is typical 



of Seyfert galaxies (eg. Bicknell et al. 1998). 
The total pressure can be expressed as: 

Ptot = Ppartides + B 2 /8lT (8) 

where 

Ppartides = Pe + Pth = (1 + fc) P e (9) 

is the total pressure in relativistic electrons P e plus the pressure in thermal particles Pth- 

The minimum pressure magnetic field is obtained by minimising the total pressure subject to 
the constraint of the observed radio surface brightness. Normally when estimating the pressure in 
the plasma of radio galaxies via minimum pressure arguments, one takes k ~ 1. For NGC 2992, 
using k = 1, an assumed spectral index of 0.7 and an electron index a = 2.4 we obtain Ppartides = 



1.1 x 10~ n dyn cm -2 , and a minimum pressure magnetic field of B = 1.5 x 10~ 5 G. Colbert et 
|aL| (1998) infer the pressure of the X-ray gas to be 1.9 x 10~ n dyn cm~ 2 for T=0.1 keV, and 
7.9 x 10~ n dyn cm~ 2 for T=l keV. This is larger than the value for Ppartides obtained in the 
minimum pressure calculation with k = 1, although only marginally so for T=0.1 keV. Therefore 
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it is possible that the conventional estimate of the minimum pressure associated with the non- 
thermal plasma is an underestimate. Moreover, we note that in order to drive a shock with velocity 
~ 350 km s _1 into gas with a pre-shock number density of ~ 1 cm -3 , we require a pressure of 
~ 1CP 9 dyn cm~ 2 . This is well above the upper X-ray estimate, although the latter is based upon 
the total extended X-ray luminosity and is therefore an average. Nevertheless the point remains 
that the pressure of the X-ray emitting gas is possibly well above the conventionally calculated 
non-thermal value. 

Therefore we consider values of the parameter k S> 1. We find that the thermal pressure, 
j^Ppartides, is equal to upper value for the pressure in the X-ray gas for k = 75. The minimum 
pressure magnetic field in this case is B = 4.1 X 1CP 5 G. If the lower X-ray pressure is relevant, we 
obtain values of k = 6 and B = 2.2 x 10~ 5 G; that are close to the minimum pressure values for 
the k = 1 case. 

Using these values of k we now estimate f e : 



fe = = [1 + (Pth/ZPeT 1 = [1 + k/2}- 1 (10) 

where e e and eth are the energy densities of the relativistic and thermal components. This 
gives / e =0.026 for k = 75, and / e =0.25 for k = 6. These low values of / e , in particular / e =0.026, 
indicate that the relativistic gas which is responsible for the radio emission may constitute only a 
small fraction of the total energy. 

Estimating the age to be t = 10 6 — 10 7 yrs from the expansion time-scale, and f a d ~ 0.5, we 
obtain r = f e f ad t = 10 4 1 - 10 5 1 yrs for / e =0.026, and r = 10 5 - 1 - 10 6 ' 1 yrs for / e =0.25. Using 
these values of r and the corresponding minimum pressure magnetic fields we can read off K1.4 
from Figure 2 in |Bickncll et al.| (1998). This gives K1.4 in the range 2 x 10~ 14 — 7 x 10 -13 Hz -1 . 
These values are larger than the observed value of K1.4 described above, implying that the energy 
associated with the radio and X-ray gas is less than the amount required to power the total optical 
emission of the ENLR. Given the uncertainties in the X-ray images and the reddening correction 
to the [O III] A5007 luminosity, we may view the results presented here as accurate to an order 
of magnitude. Even so, this result represents a significant difference compared to the traditional 
value of K4.4 ~ 10~ n which is implicitly used when the radio emission is considered as a scaled 
down version of a radio galaxy lobe. 

Another possibility is that minimum pressure conditions are not appropriate. In this case, with 
k = 1, the magnetic field required for the total particle pressure to equal the X-ray gas pressure 
is B = 3.3 x 10~ 6 G for the upper value of the X-ray gas pressure, and B = 7.5 x 10 -6 G for 
the lower value of the X-ray pressure. Using the same estimates for t and f a d as above leads to 
estimates of m,4 = 5 x 10~ 15 — 6 x 10~ 14 Hz -1 for the upper value of the X-ray pressure and 
k 1a = 3 x 10" 14 - 3 x 10" 13 Hz" 1 for the lower value of the X -ray pressure. The observed value of 
K1.4 is outside of these ranges, but only marginally so for the upper value of the X-ray pressure. 
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In the context of the theoretical work of Bicknell et al\ (1998) the properties of the optical, 
radio and X-ray emitting gas in NGC 2992 indicate that the energy associated with the radio 
plasma and the X-ray emitting gas may provide a larger fraction of the energy required to power 
the ENLRthan anticipated by previous calculations for Seyferts which assume k\a = 3x 10~ n . The 
uncertainties in the above calculations could be reduced if an accurate estimate of the temperature 
of the X-ray emitting gas was available. 



8. CONCLUSIONS 

The physical properties of NGC 2992 derived from the emission line spectra show that the 
ENLR is characterized by relatively high temperature, [S II] densities up to 400 cm -3 and a fairly 
wide range of excitation. 

In comparing the spatially mapped line ratios with an extensive shock and photoionization 
model grids, we have shown that the observations generally fall within the parameter space defined 
by either class of model. However, in terms of the consistency of the model predictions from 
one diagnostic diagram to another, the shock+precursor models provide a better outcome. The 
ionization cones are characterised by a shock velocity of 350—500 km s _1 , consistent with the line 
widths found in high resolution observations. In the inter-cone region, and about the edges of 
the ionization cones, the velocities inferred are lower, ranging down to about 200 km s _1 . Single 
component photoionization models would require a range of spectral indices in the different regions, 
which is not allowed in a model where the photons are provided by a central source. In addition, 
different line ratios demand different values of the ionization parameter in the same region, which 



is again inconsistent. The two component A M /j sequence model of pinette, Wilson, fc Storchi- 



Bergman (1996) is generally better able to account for the observed line ratios than the simple single 
component photoionization models. Indeed, the fit provided by both the shock and photoionization 
hypotheses, in particular for the higher ionization lines, would be improved by more sophisticated 
multi-component models. 

We have applied the g-test for a central ionizing source and have shown that the apparent 
increase in the rate of ionizing photons, q, with distance from the nucleus is most likely due to 
overestimation of the average density with the [S II] A6717/A6731 ratio. The most reliable values 
of q determined at the inner edge of the NLR, are however below the lower limits for the hydrogen 
recombination rate Q rec suggesting that there is a source of locally produced ionizing photons which 
may be attributed to the ionizing field generated by fast shocks. 

Despite the lack of any direct morphological indication of interaction between the radio or 
X-ray plasma with the emission line gas in NGC 2992, shocks must play an important role in this 
object. The outflow kinematics attest to this, because strong shocks with Mach numbers of order 
40 simply must occur in such dynamically active gas. Furthermore, the observation of co-spatial 



X-ray and radio emitting plasma support the entrainment theory of Bicknell et al.| (1998) in which 



the Seyfert jets are dominated by the entrained thermal plasma. In this object we find that the 
energy associated with the radio and X-ray plasma may account for a large fraction of the energy 
required to power the total optical emission of the ENLR. 

While the emission line ratios alone cannot uniquely determine shocks as the dominant exci- 
tation mechanism, the observed ratios and the velocity field and the kinematical line widths are all 
broadly consistent with shocks, and with the energetics as determined from the radio and X-ray 
plasma. 

We find no spectral evidence for significant contribution to the ionizing field or outflow from a 



starburst component, and favour the interpretation of Colbert et al. (1998) that the source of the 



outflow is an initially relativistic jet originating at the AGN. 
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Table 1. Line intensities of NGC 2992 





NW 


cone 


SE 


cone 


inter-cone 


Line 


Flux 


Flux 


Flux 


Flux 


r lux 


r lux 






(r.c) 




(r.c) 




(r.c) 


(1) 


(2) 


(3) 


(4) 


(5) 


W 


(0 


[Ne V] A3426 


560 ± 12 


81 ± 17 


61 ± 20 


130 ± 40 






[0 II] A3727 


440 ± 20 


610 ± 30 


423 ± 15 


802 ± 28 


62 ± 3 


176 ± 9 


[Nelll] A3869 


95 ± 35 


125 ± 46 


100 ± 8 


180 ± 10 


50 ± 20 


130 ± 70 


H 8 + He I A3888 


22 ± 9 


29 ± 12 










[Ne III] A3967 


31 ± 11 


40 ± 15 










[S II] A4072 






9 ± 5 


14 ± 8 






H<5 


22 ± 8 


27 ± 9 


16 ± 5 


26 ± 8 


24 ± 11 


49 ± 24 


H 7 


42 ± 4 


49 ± 4 


37 ± 6 


51 ± 8 


55 ± 13 


84 ± 21 


[0 III] A4363 


15 ± 3 


17 ± 3 


16 ± 6 


21 ± 7 






He II A4686 


28 ± 4 


29 ± 4 


32 ± 7 


35 ± 7 






H/3 


100 ± 6 


100 ± 6 


100 ± 6 


100 ± 6 


100 ± 10 


100 ± 10 


[0 III] A4959 


320 ± 20 


310 ± 20 


290 ± 40 


280 ± 40 


70 ± 5 


65 ± 5 


[0 III] A5007 


950 ± 60 


920 ± 60 


910 ± 120 


840 ± 100 


220 ± 18 


194 ± 16 


[N I] A5200 


7 ± 2 


7 ± 2 










He I A5876 


31 ± 12 


26 ± 10 


63 ± 15 


42 ± 10 


52 ± 17 


27 ± 9 


[0 I] A6300 


40 ± 3 


32 ± 3 


90 ± 14 


54 ± 8 


56 ± 11 


24 ± 5 


[0 I] A6363 


19 ± 3 


14 ± 2 


37 ± 17 


20 ± 10 






[N II] A6548 


63 ± 1 


48 ± 1 


143 ± 4 


79 ± 2 


163 ± 2 


63 ± 1 


Ha 


410 ± 5 


310 ± 3 


560 ± 10 


310 ± 8 


810 ± 10 


310 ± 4 


[N II] A6583 


189 ± 4 


142 ± 3 


430 ± 10 


237 ± 7 


490 ± 7 


185 ± 3 


[S II] A6716 


123 ± 6 


91 ± 5 


230 ± 20 


120 ± 10 


220 ± 32 


79 ± 11 


[S II] A6731 


95 ± 5 


71 ± 4 


180 ± 20 


98 ± 9 


174 ± 25 


62 ± 9 


[S III] A9069 


67 ± 5 


38 ± 3 


83 ± 42 


25 ± 13 


251 ± 27 


37 ± 4 


[S III] A9532 


197 ± 11 


109 ± 6 


300 ± 70 


90 ± 20 


512 ± 44 


68 ± 6 



Note. — All intensities are given relative to H/3=100. (r.c) indicates reddening corrected 
intensities. 
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Table 2. 



Region 


L 40 ([O III]) 


L 40 (Ha) 


logloQrec 


NW cone 


5.01 


1.9 


53.2 


SE cone 


5.63 


3.7 


53.5 


NW cone (red corrected) 


14.4 


4.2 


53.5 


SE cone (red corrected) 


76.9 


29.9 


54.4 


Entire EELR 


11.7 


6.70 


53.0 


Entire EELR (red corrected) 


96.0 


40.1 


53.8 



Note. - Luminosities of the NW, SE and inter-cone regions of 
NGC 2992 in units of 10 40 ergs s _1 , and the inferred rate of hydrogen 
recombinations Q rec (s -1 )- 
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Fig. 1— Optical and radio images of NGC 2992: a) [0 III] A5007 narrow band image, b) Ha+[N II] 
narrow band image, c) blue continuum image, d) blue continuum image overlaid with the 20cm 
radio map from [Hummel et aL] (1983) (contours: -2, 2, 4, 7, 10, 20, 30, 50, 70, 90, 110 mJy per 
beam; beam=6") e) [O III] narrow band image overlaid with 20cm radio contours, f) [O III] overlaid 
with the 6cm radio map from |Ulvestad and Wilson| (1984) (contours: 3, 6, 9, 15, 30, 50, 70, 



of the peak of 7.2 mJy per beam; beam=0.5"). 

Fig. 2. — Slit positions used for the long-slit observations of NGC 2992 overlaid on the [O III] A5007 
contour map with contour levels shown at 0.5, 1, 2, 5, 10, 20, 40, 60, 80, 100% of the peak value. 

Fig. 3. — a) Average spectrum of the NW cone of NGC 2992, b) Average spectrum of the SE cone 
of NGC 2992, c) Spectrum of the inter-cone region of NGC 2992. 

Fig. 4. — Blue spectra of the NW and SE cones of NGC 2992. The top curve shows the observed 
spectrum overlaid with the template used to subtract the underlying stellar contribution. The 
lower curve shows the template subtracted spectrum. 

Fig. 5. — Spatial variation of line ratios in NGC 2992. The open circles indicate line ratios from 
positions in the NW cone, solid squares indicate points in the SE cone, and open triangles indicate 
the inter-cone region as shown in Figure 6. 

Fig. 6. — Spatial locations of the extracted spectra overplotted on the [O III] A5007 contour map. 
The contour levels are as in Figure 2. The open circles indicate positions in the NW cone, solid 
squares indicate points in the SE cone, and open triangles indicate the inter-cone region. 

Fig. 7. — Ha + [N II] A6548,6583 profiles in the spectra of the NW cone (top), SE cone (center), and 
the inter-cone regions (lower). In each panel we show the total fit overlaid on the data template 
subtracted data. The lighter weight curves show the fitted galaxy template and the individual 
Gaussian component fits to the Ha and [N II] emission lines. 

Fig. 8. — Electron densities n e inferred from the [S II] A6717/A6731 ratio in spectra extracted from 
various positions in the ENLR of NGC 2992. The plotting symbols represent the different regions 
of the ENLR as shown in Figure 6. 

Fig. 9. — Ionization parameter estimated from the [O IIl]/[0 II] line ratio at each point in the 
spectral map of NGC 2992. The plotting symbols represent the different regions of the ENLR as 
shown in Figure 6. 

Fig. 10. — [O III] A5007/Ha versus [N II] A6583/Ha. The observed line ratios are shown for each 
point in the spectral map of the ENLR of NGC 2992. The plotting symbols represent the different 
regions of the ENLR as shown in Figure 6; open circles denote points from the NW cone, solid 
squares denote points from the SE cone, and open triangles indicate the inter-cone region. In 
panel (a) we over-plot the shock, and shock+precursor model grids and also include the average 
Ict error bar and a vector which indicates the direction of reddening correction. The shock and 
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shock+precursor model grids are labelled with shock velocity in units of km s , and with the 
magnetic parameter B/n 1 / 2 = 0—4 in units of /j,G cm 3 / 2 . Panel (b) shows the single component 
photoionization (a=-l. 0,-1.4 and -2.0) [/-sequences plotted for n = 100 cm~ 3 and n = 1000 cm" 3 
with grid lines of constant U every 0.25 dex, and labelled with logio^. The Am/i sequence is also 
shown and is labelled with A M /i an d has tick marks at 0.2 dex intervals. 

Fig. 11. — [O III] A5007/Ha versus [N II] A6583/Hq. The data points and model grids are as 
described for Figure 10. 

Fig. 12.— [O III] A5007/H/3 versus [S II] AA6717,6731/Ha. The data points and model grids are 
as described for Figure 10. 

Fig. 13.— [O III] A5007/H/3 versus [O II] A3727/[0 III] A5007. The data points and model grids 
are as described for Figure 10. 

Fig. 14.— [N II] A6583/Hq versus [S II] AA6717,6731/Ha. The data points and model grids are as 
described for Figure 10. 

Fig. 15.— [O I] A6300/[O III] A5007 versus [O II] A3727/[0 III] A5007. The data points and model 
grids are as described for Figure 10. 

Fig. 16.— [S III] AA9069,9532/Ha versus [S II] AA6717,6731/Hq. In panel (a) the shock-only grid 
is shown with the dashed line. The data points and model grids are otherwise as described for 
Figure 10. 

Fig. 17. — [S II] AA6717,6731/Ha versus He II A4686/H/3. In panel (a) the shock-only grid is shown 
with the dashed line. The line ratios for the average spectra of the NW and SE cones are shown as 
the open circle, and the solid square respectively. The model grids are as described for Figure 10. 

Fig. 18.— He II A4686/H/3 versus [O II] A3727/[0 III] A5007. The line ratios for the average spectra 
of the NW and SE cones are shown as the open circle, and the solid square respectively. The shock, 
shock+precursor, [/-sequence and A M /i sequence models are shown and are labelled as described 
in Figure 10. 

Fig. 19.— [Ne V] A3426/[Ne III] A3869 versus [O III] A5007/H/3 ratio. The line ratios for the 
average spectra of the NW and SE cones are shown with the open circle, and the solid square 
respectively. The shock, shock+precursor, [/-sequence and A M /j sequence models are shown and 
are as described in Figure 10. 

Fig. 20. — Temperature sensitive diagnostic diagram: [O III] A4363/[0 III] A5007 versus He II 
A4686/H/3. The Rom temperature scale is shown on the left side of the diagram. The line ratios 
for the average spectra of the NW and SE cones are shown with the open circle, and the solid square 
respectively. The shock, shock+precursor, [/-sequence and A M /i sequence models are shown and 
are as described in Figure 10. 
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Fig. 21. — Densities (a), and Ionization parameters (b) for points within the NW and SE cone 
regions that were used in the q-test calculations. 

Fig. 22. — Q versus radius, and n e: Q versus radius for three different de-projections of the points 
within the ionization cones of NGC 2992. The diagrams in the left panels show Q versus radius 
for de-projection onto the near side for both cones (top), de-projection onto the far side of the 
cones (middle) and with no de-projection (bottom). The diagrams in the right panels show the 
corresponding plots of n e< Q versus radius. The plotting symbols represent the different regions of 
the ENLR as shown in Figure 6. The lower useful limit for the [S II] density diagnostic is shown 
with the dashed line in the diagrams in the right panels. 

Fig. 23. — a) Outflow velocity field overplotted on the [O III] A5007 contour map. The size of 
the circle at each point indicates the difference in velocity of the wind component relative to the 
rotational component. The key in the top right hand corner indicates the magnitude of this velocity 
in km s^ 1 . Solid points correspond to velocities that are blue shifted relative to the rotation, and 
open circles indicate red shifted velocities, b) FWHM map of the major outflow component. The 
image contours are the same as in Figure 2. 

Fig. 24. — Excitation versus line width diagrams, (a): [O III] A5007/H/3 versus the FWHM of 
the kinematical components. The kinematical component identified with the galaxy rotation is 
plotted as the solid circles, the main outflow components are plotted as open squares and the 
extra component (detected at a few locations) is plotted as open triangles, (b): [N II] A6583/Ha 
versus FWHM. As in (a), the square points indicate the main outflow component, and the triangles 
indicate the extra component. In (b) we make the distinction between the cones by using filled 
symbols for the SE cone and open symbols for the NW cone. In both (a) and (b) the solid curve is 
the shock -l-precursor model prediction where the FWHM is directly related to the shock velocity, 
and the dashed curve assumes that the FWHM is associated with the projection of randomly 
oriented shocks so that V s h oc k ~FWHM/\/3. The uncertainty in the FWHM is generally less than 
100 km s _1 . 
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